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I. INTRODUCTION 



The study of low-energy QCD provides insight into the non-perturbative dynamics of 
quarks and gluons. One of the most fundamental observables is the spectrum of the theory. 
Quarks and gluons are confined into color neutral hadronic states that are seen in nature. 
Understanding the mass spectrum of these lowest-lying states directly from QCD remains 
an on-going challenge. 

A model independent tool to study low-energy QCD is through chiral perturbation theory 
(xPT). This is an effective theory written in terms of low-energy degrees of freedom, e.g. the 
octet mesons in SU(3) flavor are assumed to be the pseudo-Goldstone bosons that appear 
from chiral symmetry breaking. In this effective theory observable quantities receive both 
long-range and short-range contributions. The long-range contributions arise from non- 
analytic meson loops, while the short-range physics is encoded in a number of low-energy 
constants (LECs). Symmetry constrains the number of such constants but their values must 
be determined from experiment or from theoretical calculations. 

Lattice QCD can provide first principles numerical determination of QCD observables, in 
particular the mass spectrum of baryons. There have been numerous calculations of baryon 
masses m quenched QCD (QQCD) and a few in QCD 0, Q 3 3 Partially quenched 
QCD (PQQCD) calculations in the baryonic sector are so far very limited. A problem that 
haunts lattice calculations is that computing power currently restricts the simulations to 
quark masses that are larger than the physical light quark masses. Therefore, to make 
physical predictions it is necessary to extrapolate from the heavier lattice quark masses to 
those found in nature. For QQCD, where the fermionic determinant that arises from the 
path integral is set to a constant, quenched chiral perturbation theory (Q^PT) has been 
developed to aid in the extrapolation p, 

HiSBQQ- There is no general connection 
of QQCD observables to QCD because QQCD does not have an axial anomaly. This feature 
of the quenched approximation leads to new operators in QxPT that have no analogue in 
%PT. Moreover the LECs of QxPT are numerically different than in xPT. 

The problems of the quenched approximation can be remedied by using partially quenched 
QCD (PQQCD). In PQQCD contributions from sea quarks are retained and the fermionic 
determinant is hence no longer a constant. Additionally the masses of sea quarks are varied 
independently of the valence quarks. By efficaciously giving the sea quarks larger masses it 
is much less costly to calculate observables than in ordinary QCD. The low-energy effective 
theory of PQQCD is partially quenched *PT (PQx?T) El H H H H B i 



Since PQQCD retains an axial anomaly, the singlet field can be integrated out. Therefore 
the LECs of xPT appear in PQxPT. By fitting PQ^PT to PQQCD lattice data one can 
determine the LECs and thereby make physical predictions fo r Q CD. There has been much 
activity recently in calculating baryon properties in PQ^PT |22l I23I Q 0, 27, H^, 
3oL 31, 3^, 33|. There has also been interest in using chiral effective theories to ex trap olate 
quenched and unquenched lattice data on the baryon masses, see for example jH 35[ 36| . 

In this work we calculate the masses of the decuplet baryons to next-to-next-to leading 
order (NNLO) in xPT. We also provide the NNLO calculation of the decuplet baryon masses 
in PQ^PT. These calculations are performed in the isospin limit of SU(3) and in the partially 
quenched calculation we use three valence, three ghost and three sea quarks, with two of the 
sea quarks degenerate. The expressions we derive can be used to extrapolate lattice QCD 
and partially quenched lattice QCD data to the physical quark masses. Additionally the 
LECs can be determined up to 0(m^), where m q represents a quark mass. 



2 



This paper has the following organization. First in Section |HJ we review the basics of 
%PT including the inclusion of the octet and decuplet baryons to leading and next-to-leading 
order in the heavy baryon expansion. In this work we use rank-three flavor tensors for both 
the octet and decuplet baryons. Next in Section lllll we calculate the decuplet masses to 
next-to-next-to-leading order in xPT. We work in the SU(3) flavor group and in the isospin 
limit. In Section HVI we review PQ^PT in a way that parallels our discussion of xPT, then 
in Section we calculate the decuplet masses at next-to-next-to-leading order in PQ^PT. 
Finally a summary (Section IVIj) ends the paper. 



II. HEAVY BARYON CHIRAL PERTURBATION THEORY 

In this section, we briefly review chiral perturbation theory. We start first in the meson 
sector and then review heavy baryon chiral perturbation theory to leading and next-to- 
leading order in the inverse baryon mass. For the baryons, we differ from the standard 
formulation by embedding the octet states into rank-three flavor tensors. 



A. Pseudo-Goldstone bosons 



For massless quarks, the QCD Lagrangian exhibits a chiral symmetry, SU (3)l®SU(3)r® 
U(l)v, which is spontaneously broken down to S77(3)y <8>{7(l)y. Chiral perturbation theory, 
the low-energy effective theory of QCD, emerges by expanding about the physical vacuum 
state. Without the explicit breaking of chiral symmetry provided by the quark mass term 
of the Lagrangian, eight particles (the pions, kaons and eta) would emerge as the Goldstone 
bosons of the broken SU {3) a symmetry. Given that the quark masses are small compared 
to the scale of chiral symmetry breaking, however, the lowest lying mesons emerge as an 
octet of pseudo-Goldstone bosons. 

The pseudo-Goldstone bosons are collected in an exponential matrix 



exp 



2 i $ 





(>°+ 




7T + 


$ = 


7T 










K° 



K + \ 

-iky 



With the above convention, the pion decay constant / is 132 MeV. The effective Lagrangian 
describing the dynamics of these mesons at leading order in xPT is |37( 



£ = L tr (dPYtdfX) + A tr (m ] Q Y, + m Q S 1 



(2) 



Above, leading order in the power counting is 0{mq) and hence Q{mq) ~ 0(q 2 ), where q 
is an external pion momentum. In the isospin limit, the quark mass matrix tuq is given by 



m Q = diag(m n , m u , m s 



(3) 



Expanding the Lagrangian Eq. (J2J to leading order, one finds that mesons with quark content 
QQ' are canonically normalized when the meson masses are given by 



m 



QQ 1 



4A 
7 



2 ( m Q + m Q' 



(4) 



3 



B. Baryons 



Besides the mesons, there are also three quark states in the spectrum of QCD. To include 
these octet and decuplet bar yon s systematically into chiral perturbation theory, we use heavy 
baryon xPT (HBxPT) H H ffl S3 . For the octet, baryon fields B (x) are redefined in 
terms of velocity dependent fields B v (x), 

Bv{x) = ±±le lM * v - x B(x) , (5) 

where is the four-velocity of the baryon, B. This field redefinition corresponds to param- 
eterizing the baryon momentum as 

p M = M B v^ + k^ (6) 

where k^ is the residual momentum. This efficacious redefinition eliminates the Dirac mass 
term for baryons 



B (i? -M B )B = B v i$B v + O . (7) 

From Eq. J5J), it is easy to verify that derivatives acting on B v bring down powers of the 
residual momentum k. Thus, higher dimension operators of the heavy baryon field B v are 
suppressed by powers of Mb and a consistent derivative expansion emerges. Henceforth we 
shall omit the subscript v from all baryon fields treating them implicitly as heavy baryons. 

In this work we embed all baryons in rank-3 flavor tensors. The convenience of this 
choice will become readily apparent when we generalize to PQ%PT. The SU(3) matrix of 
the lowest-lying spin-| baryon fields is 



p 



B=| ^A-^£° n \- (8) 

These baryon states are embedded in the tensor B^ in the following way [1 

B ijk = {^ijl^k + 6 ikl~Bj) ■ (9) 

The flavor tensor has the symmetry properties 

B^k = Bikj and B^ + Bjik + Bkji = 0. (10) 

When the spin-| decuplet baryons T are included in the theory, an additional mass 
parameter A appears. This parameter is the leading-order mass splitting between the octet 
and the decuplet in the chiral limit and must be included in the power counting. We 
treat A as 0(q), where q is a typical small pion momentum. The spin-| decuplet baryons 
can be described by a Rarita-Schwinger field, (T^)^, which is totally symmetric under the 
interchange of flavor indices. We employ the normalization convention in which T 111 = A ++ . 
The heavy baryon Rarita-Schwinger field satisfies the constraints, v • T — S ■ T — 0, where 
5^ is the covariant spin-vector. 
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The Lagrangian to leading order in the 1/Mb expansion can be written as 

£ = (BivDB) +2a M (BBM+) + 2{3 M (BM+B) + 2a M (BB) tv(M+) 
- (T" [iv ■ D - A] T„) + 2 lM (T^M+T^) - 2o M tr(M+) 
+2a (BS^BA^ + 2f3(BS^A ll B) + 2H (T u S^A^) 

+ [(T%B) + (SA„T")] . (11) 

Above, _D M is the chiral-covariant derivative which acts on the B and T fields as 

(D»B) ijk = d»B ijk + (V^ l B ljk + (Vn'Buk + (W) k l B ijh (12) 
The vector and axial-vector meson fields appearing in the Lagrangian are given by 

V^^td^ + gd^), A, = ^8^-^0,0, (13) 

and 

m+ = \ (e w f + ■ ( 14 ) 

In Eq. (111)1 the brackets ( ) denote a contraction of the flavor indices and are defined 
in (l2|. Such contractions ensure the proper transformations of the field bilinears under 
chiral transformations. To compare with the coefficients used in the standard two-index 
baryon formulation |38l 3^, 4(], 41 1, it is straightforward to show 



and 



a = l D + 2 F, (3 = -^ D + F, (15) 

2 5 
olm = ~bD + 2&f , Pm = - -bD + b F , <Jm = b D - b F + a . (16) 



C. Higher dimensional operators 

The Lagrangian in Eq. (jllj) contains some, but not all, terms of 0(q 2 ). To calculate the 
decuplet masses to 0(71%) we must include all 0(q 2 ) relevant operators, which contribute 
to the mass via loops, and all 0(q 4 ) relevant operators which contribute at tree level. The 
Lagrangian also includes operators of 0(q 3 ), but they do not contribute to the self-energy. 

The baryon mass is treated as Mb ~ A x . As the LECs are a priori unknown, we can 
combine the 1/Mb and 1/A X expansions into one expansion in powers of 1/A X . There is one 
exception: constraints from reparameterization invariance (RI) determine the coefficients of 
some of the higher dimension operators arising in the 1/Mb expansion |!2,|43||. Thus these 
1/Mb corrections must be kept distinct to insure the Lorentz invariance of the heavy baryon 
theory to a given order. 

The heavy baryon momentum parameterization in terms of ^ and k 11 in Eq. (jUJ) is unique 
only up to 1/Mb corrections. When the velocity and residual momentum are simultaneously 
transformed in the following way 

V ^ V + TT , k^k-e, (17) 
Mb 
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the momentum in Eq. © is unchanged. Reparameterization invariance requires the 
effective Lagrangian to be invariant under such transformations and thus ensures the theory 
is Lorentz invariant to a given order in 1/Mb- Furthermore, utilizing RI has non-trivial 
consequences as it connects operators of different orders in the 1/Mb expansion and thereby 
exactly fixes the coefficients of some of the higher dimensional operators with respect to the 
lower ones. We find that the fixed coefficient Lagrangian is 



£= - 



- D* 
B — ^B 

2M B 



_ n 2 

T»—±-T„ 



2M. 



B 



n 



T M — — — v ■ AT,, 



M 



S ■ %D 
T M v ■ A——T„ 



B 



M. 



B 



18) 

where D\ = D 2 — (v ■ D) 2 and we have kept only terms relevant to the calculation of decuplet 
masses. 

The Lagrangian contains additional 0(q 2 ) operators as well as 0(q 4 ) operators which 
are invariant under the £77(3) chiral transformations. In part, these operators absorb some 
effects of the unphysical, off-shell degrees of freedom 44]| . The operators relevant to the self 
energy of the decuplet baryons are 



t A Tf (AJTtfTfa + t*Tf(A u )/(An/T{: rk + t A pyT") tr(A u A») 



4tt/ 



+t A Tf{A»A v )/T^ k + t§Tf W{A„)/T», jlk + t§ (T,T») tr(AM„) 

+tl A Tf(v ■ A v ■ A)/TP jk + tl A Tf{v ■ A)/(v ■ A)/T$ jth + (T M T^) tr(v ■ Av ■ A) 

+t A 1 I Tf(M + M + )/T» + t^Tf(M + )/(M + )/T^ k + tf {T»T») tv{M + M + ) 



+tf (T^M + T^)tT(M+)+t^ (T^)tv(M+)tT(M. 



(19) 



All the LECs tf, tf, , and tf' 1 are dimensionless. In principle, additional 1/Mb operators 
with the same chiral symmetry properties as those contained in Eq. can be generated. 
However, these 1/Mb operators do not have their coefficients constrained by RI, therefore 
they shall be absorbed into the definition of the various tf, tf, t^ A , and tf 1 . Since the flavor, 
spin and Lorentz structure of these omitted 1/Mb operators is identical to those above, we 
are guaranteed that the values absorbed in the LECs of Eq. (fT9"j) remain the same for all 
processes and thus the determination of the above LECs is in that sense universal. 

Including the decuplet fields in x?T requires additional operators involving A/A x because 
A is a chiral singlet. Thus arbitrary functions f(A/A x ) can multiply any term in the 
Lagrangian without changing the properties under chiral transformations. To the order 
we are working, all constants in the calculation must be arbitrary polynomial functions of 
A/A x , and expanded out to the appropriate order. For example 



7m — > 1m 



A 



7m 



1 + 7i 



A 

AT 



A 2 



A3 



(20) 



Furthermore at this order the decuplet-octet mass splitting in the chiral limit is a polynomial 
function of A/A x . We shall not explicitly write down the operators that contribute to the A 
dependence because determination of their LECs requires the ability to vary A. Additionally, 
the LECs are also implicit functions of fi, to absorb the divergences arising from the loop 
integrals. 
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FIG. 1: One loop graphs which give contributions to Mj,. ■ The single, double and dashed lines 
correspond to octet baryons, decuplet baryons and mesons, respectively. The filled squares denote 
the axial coupling given in Eq. 



TABLE I: The tree level coefficients m?, (m 2 )T, and (mm')? in xPT and PQxPT for decuplet 
states T. 





mj- 


(m 2 ) T 


(mm')x 


A 


3m u 


3m 2 u 


^ 3m2 




2m u + m s 


2ml + m s 


ml + 2m u m s 




m u + 2m s 


ml + 2m 2 s 


2m u m s + m 2 




3m s 


3m 2 s 


3m 2 



III. DECUPLET MASSES IN xPT 



The masses of the octet and decuplet baryons have been investigated considerably in 
XPT 0, H H H H 5<j . Here we calculate the masses of the decuplet baryons to 
NNLO in xPT. The mass of the i th decuplet baryon in the chiral expansion can be written 

as 

M Ti = Mo (//) + 4 J (//) + Mg /2) (,x) + jwg> (/i) + . . . (21) 

Here, Mq (//) is the renormalized mass of the decuplet baryons in the chiral limit which is 
independent of tuq and also of the T,. is the contribution to the i th decuplet baryon of 
the order mg', and /i is the renormalization scale. For this calculation we use dimensional 
regularization with a modified minimal subtraction (MS) scheme. 

In calculating the masses, the leading dependence upon itlq arises from the terms the 
Lagrangian Eq. (fTTj) with coefficients and Wm- The 0(m i Q 2 ) contributions arise from the 
one-loop diagrams shown in Fig. ^ which are formed from the operators in the Lagrangian 
with coefficients 7i and C. The 0(m 2 Q) contributions arise from the one-loop diagrams shown 
in Fig. |21 from the tree-level contributions of the operators with coefficients, tf f , and from 
the NLO wave-function corrections. 

We find that the leading-order contributions to the decuplet masses are 

2 

= - 7m m T - 2a M tr(m Q ), (22) 



where the coefficients ttit appear in Table UJ The next-to-leading order contributions are 

r (3/2) _ 5ft 2 ^ T o 
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TABLE II: The coefficients AT and BT in xPT f° r decuplet states T. 



A 



TT 

~5~ 

6 
4 
9 
1 
6 





.4 



T 



if 

1 

3 





1 

6 
2 
3 



7T 

~T~ 

3 
5 

1 
3 





B 



■T 

K 

2 
3 
4 
9 
2 
3 
4 
3 



1 

3 
1 

3 





\ 

=□= 



(2) 

FIG. 2: One-loop graphs which give contributions to Mj, . Single, double and dashed lines 
correspond to octet baryons, decuplet baryons and mesons, respectively. The filled squares and 
filled circles denote axial couplings and insertions of A4 + given in Eq. The empty squares 

denote insertions of fixed 1/Mb operators from Eq. ()18|) . and the empty circles denote two-baryon- 
axial couplings defined in Eq. ()19|) . 



where the function T is denned by 
jF(m, 5, /i) 



V S 2 — m? log 



5 m 2 log 

2 6 



/// 



5 — \J5 2 — m 2 + is 
5 + y/5 2 — m 2 + is 



8 log 



(24) 



In Eq. (|25j) the sum on runs over loop mesons with the mass m^. The coefficients and 
above are the sums of squares of Clebsch-Gordon coefficients. The numerical values of 
these coefficients depend on the decuplet state T and are listed in Table |Il] for loop mesons 
with mass m^. 

Finally the next-to-next-to-leading contributions to the decuplet mass are 
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TABLE III: The coefficients and C$ in xPT. Coefficients are listed for the decuplet states T 
and grouped into contributions from loop mesons with mass m^. is identical for each decuplet 
state. 





7T 


K 


V 


A 

S* 

n- 


3m u 
2m u 
m u 



m u + m s 
|(m„ + m s ) 
\{m u + m s ) 
2(m u + m s ) 


l(m u + 2m s ) 
\{m u + 8m s ) 

4 




6m u 


\(m u + m s ) 


|(m„ + 2m s ) 



TABLE IV: The coefficients -D^ and in xPT. Coefficients are listed for the decuplet states T 
and grouped into contributions from loop mesons with mass m^. is identical for each decuplet 
state. 





TT 


K 


7] 


A 


3 
2 


1 


1 

6 


S* 


1 


4 
3 


1 

3 




1 


5 


1 




2 


3 


2 


n- 





2 


2 
3 




3 


4 


1 



Mr 



(2) 



(Z - 1)M 

i r i 



(i) 



2 7 M 



l -t™ (mm') T + ff tr(mj) + m T tr(m Q ) + [tr(m Q )] 2 | 



(4tt/) 



2 E + (4^2 E C * £ ( m ^' ^ 



(4tt/) 



4(4tt/) 
5 7i 2 



(4tt/) 2 M b 



£4 



15 C 2 



\m 2 o A Mm Q ) 

3(4tt/)2 
10H 2 7 m 



£4 



'/./ \ 26 2 
£(m^,/i) + — m 



16 (4vr/) 2 M B 
3C 2 cr M tr(m Q ) 



Cim^y) - — mj 



9(4tt/) s 



E F J 



'/./ \ 26 2 



3C 2 



2(4tt/) : 



E^ 

(25) 



(47T/) 2 

^Gj^m^-A,//). 



Here the tree-level contributions are expressed in terms of coefficients (m 2 )^, and 
(mm')r which are listed in Table |U Above the wavefunction renormalization Z is given 
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TABLE V: The coefficients in ^PT. Coefficients are listed for the decuplet states T and grouped 
into contributions from loop mesons with mass tha,. 





7T 


K 


V 


A 

S* 
"* 


1 

2 
1 

6 






2 
3 
2 
3 


1 

6 

1 
6 




n- 








2 
3 



TABLE VI: The coefficients in xPT- Coefficients are listed for the decuplet states T and 
grouped into contributions from loop mesons with mass m^. 





7T 


if 


rj 


A 


5 

2 WI"U 


|(2m u + m s ) 




£* 


§(2m u + m s ) 


|(2m u + m s ) 







\{m u + 2m s ) 


i(4m u + 5m s ) 


g(m u + 2m s ) 


Q~ 





|(m u + 2m s ) 


2m s 



by 



5H 2 



3(4tt/) : 



£4 



£(?7l0,/i) + — /»" 



26 

7 

15 



3C 5 



2(4tt/)' 



J2^J(m^-A,fi). (26) 



The above equations (j2*B]) and ()26|) also employ abbreviations for non-analytic functions 
arising from loop contributions. These functions are 



C(m, fi) = m 2 log 
£(m, jj) = m log 



m 



9 

m 
V 2 



J(m, 5, fj) = (m 2 - 2S 2 ) log + 28V S 2 - m 2 log 



5 — V 8 2 — m 2 + 



(27) 
(28) 
(29) 



As in Eq. (|23jl . the sums over in Eqs. (125)1 and (j2Bj) run through loop mesons with mass 
m$. The coefficients in these sums are themselves the sums of squares of Clebsch- Gordon 
coefficients and/or quark mass insertions. These are listed in Tables ITTTMVIH To be clear, 
the above expressions Eqs. (|2*2*|) . (|2~3|) . and (|2*3)l are quark mass expressions, i.e. , the meson 
masses in these expressions are merely replacements for the quark masses via Eq. (j3J). To 
utilize these expressions for chiral extrapolations of lattice QCD data, one must write the 
quark masses in terms of physical meson masses (as determined on the lattice). To work 
consistently to 0(m 2 ), one needs to use the NLO expressions for the meson masses in the 
LO results in Eq. <^ 
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TABLE VII: The coefficients in xPT. Coefficients are listed for the decuplet states T and 
grouped into contributions from loop mesons with mass m^. 





7T 


K 


V 


A 


\m u (a M + Pm) 


l[m u (5a M + Wm) 
+m s (a M + 4/3m)] 





S* 


±[m u (19a M + 22(3 M ) 


^[m u (7a M + 10/3m) 


i[m u (5a M + 2/3m) 




+m s (lla M + 80m)] 


+m s (5a M + 2/?m)] 


+m s (a M + 4/3m)] 


"* 


\[m u {a M + 4/3m) 


|[2m u (a M +/3m) 


|K(aM + 4/?m) 




+m s (5a M + Wm)] 


+m s (aM + /?m)] 


+m s (5aM + 2/?m)] 


n~ 





|[m u (a M + 4/3m) 
+m s (5a M + 2/?m)] 






IV. PQxPT 



In PQQCD, the quark part of the Lagrangian is 

9 

£= ^g^-m Q )/Q fe . (30) 
j,k=i 

This differs from the SU (3) Lagrangian of QCD by the inclusion of six extra quarks; three 
bosonic ghost quarks, (u, d, s), and three fermionic sea quarks, (j, /, r), in addition to the light 
physical quarks (u, d, s). The nine quark fields transform in the fundamental representation 
of the graded SU(Q\3) group. They have been accommodated in the nine-component vector 

Q j = (u,d,s,j,l,r,u,d,s). (31) 
The quark fields obey the graded equal-time commutation relation 

QfMQf (y) - (-i)^gf (y)Q?(x) = rt^ 3 (x - y), (32) 

where a, (3 are spin and i,j are flavor indices. Analogous graded equal-time commutation 
relations can be written for two Q's and two Q^s. The grading factors 

fl for h = 1,2,3, 4,5, 6 
lh \ for k = 7, 8, 9 1 J 

take into account the different fermionic and bosonic statistics of the quark fields. In the 
isospin limit the quark mass matrix of SU(6\3) is given by 

m Q = diag(m n , m u , m s , rrij, mj,m r , m u , m u , m s ). (34) 

Because the ghost quark masses are identically equal to the valence quark masses there is 
an exact cancellation in the path integral between the valence quark determinant and the 
ghost quark determinant. The sea quark determinant is unaffected. Thus in PQQCD, one 
has the ability to vary the valence and sea quark masses independently. QCD is recovered 
in the limit rrij — > m u and m T — > m s . 
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A. Pseudo-Goldstone Mesons 



For massless quarks, the theory corresponding to the Lagrangian in Eq. (JSUj) has a graded 
SU(6\3)l<S> SU(6\3)n<%)U(l)y symmetry which is assumed to be spontaneously broken down 
to SU(6\3)v <8> U(l)v in analogy with QCD. The effective low-energy theory obtained by 
perturbing about the physical vacuum state of PQQCD is PQxPT. The result is 80 pseudo- 
Goldstone mesons with dynamics described at leading order in the chiral expansion by the 
Lagrangian 

f 2 

£ = ^-str (d^E^E) + A str (m ? E f + m+E) + a^%d^% - m 2 $ 2 , (35) 



where 



and the meson fields appear in 



ew[ 2 -f)=e, (36) 



The operation str() in Eq. (j33|) is the supertrace over flavor indices. The quantities a$ and 
mo are non- vanishing in the chiral limit. M and M are matrices containing bosonic mesons 
(with quantum numbers of qq pairs and qq pairs, respectively), while \ an d X* are matrices 
containing fermionic mesons (with quantum numbers of qq pairs and qq pairs, respectively). 
The upper 3x3 block of M is the usual octet of pseudo-scalar mesons and the remaining 
components are mesons formed with one or two sea quarks, see e.g. |22j |. 

The flavor singlet field is defined to be $o = str($)/v / 6- PQQCD has a strong axial 
anomaly U (1)a and therefore the mass of the singlet field mo can be taken to be the order of 
the chiral symmetry breaking scale, mo — > A x Jljj. In this limit, the rj two-point correlation 
functions deviate from their form in x?T. For a,b — u, d, s, j, I, r, u, d, s, the r] a r]b propagator 
with 2 + 1 sea-quarks at leading order is 



2 



VaVb ~ q 2 - m 2 Va +ie 3 (g 2 - m 2 a + it) (q 2 - m 2 b + ie) (q 2 -m 2 x +ie)' ' 



where 



e a = (-l) 1+ ^. (39) 

The mass m xy is the mass of a meson composed of (anti)-quarks of flavor x and y, while the 
mass mx is defined as m 2 x = l(m 2 j + 2m 2 rr ). The singlet propagator can be conveniently 
rewritten as 

GrtaVb = e a$abPa + ^bHab (Pa, A, Px) , (40) 
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where 



q 2 — m 2 + it 



Ph 



q 2 — m 2 b + ie 



P 



x 



q 2 — m x + ie 



n ab (A,B,C) 



1 

3 



(m% - m 2 a ) (m 2 rr - m 2 a ) ^ _ (m% - m 2 J (m 2 rr - m 2 J 



( 



'la 



m 



) ( m Va ~ m x) ( m r, 



ihaa) (jn 2 x ~ m 



4 



x 



m 2 a ) (m 2 x 



C 



m 



>ib. 



m. 



vi, 



)( 



m 2 x ) 



B 



(41) 



B. Baryons 

In PQxPT the baryons are composed of three quarks QiQjQk where i — k can be valence, 
sea or ghost quarks. One decomposes the irreducible representations of SU (6\3)v into irre- 
ducible representations of SU(3) va i®SU(3) sea <g>SU(3)ghost®U(l). The method for including 
the octet and decuplet baryons into PQxPT is to use the interpolating field [l2[ H2] 



BJ* ~ {QT a Qf b QT - QT a QTQi b ) (42) 

Under the interchange of flavor indices, one finds [l2| 

B ijk = (-) 1+r > jm B ikj and B ijk + {-) 1+ ^B jik + (-f^m+wu+w^.. = . (43) 

We require that B^k = B^k, defined in Eq. (0), when the indices, k are restricted to 1 — 3. 
Thus the octet baryons are contained as an (8, 1, 1) of SU(3) va i®SU(3) sea ®SU(3)gf l0 st®U(l) 
in the 240 representation. In addition to the conventional octet baryons composed of valence 
quarks, Bij k also contains baryon fields composed of sea and ghost quarks. In this work we 
only need states of the 240 which contain at most one sea or ghost quark, and these states 
have been explicitly constructed in 2^ . 

Similarly, one can construct the spin-| baryons which make up the 138, and have an 
interpolating field 

t%? - (Q? a Qf b Ql c + Qi ,b Q] ,c Q k a + Q? c QT a Ql' b ) (CV% 7 . (44) 

Under the interchange of flavor indices, one finds that 

%jk = (-) 1+r >^T Jlk = {-f^% kj . (45) 

We require that %j k = T^ k) when the indices i,j, k are restricted to 1 — 3. Under SU(3) va i ® 
SU(3) sea ® SU(3) g host ® U(l) they transform as a (10, 1, 1). In addition to the conventional 
decuplet resonances composed of valence quarks, %j k contains fields with sea and ghost 
quarks. As with the 240, for our calculation the required states of the 138 have been 
constructed in 22| . 



To write down the PQ^PT Lagrangian, we must also include the appropriate grading 
factors in the contraction of flavor indices. These are included in the () notation as was 
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originally defined in 0, 22|| The leading order PQxPT Lagrangian is given by 

£ = (BivD B) + 2a M (BBM+) + 2(3 M (BM+B) + 2a M (BB) str(M- 
- (7* [ iv ■ D - A ] Tj + 2 7M (T"A4 + Tj - 2cf A/ (T%) str(A4+) 
+2a (H^B^) + 2/3 (ESM^B) + 2ft (T v S' l A l X) 



(46) 



The low-energy constants appearing above have the same numerical values as those in %PT. 

At higher orders, the situation is quite similar to the %PT case considered in Section ITT1 
Recall that at higher orders, the Lagrangian can contain arbitrary functions of A/A x . We 
take this into account by implicitly treating the leading order coefficients as functions of 
A/A x expanded out to the required order. To enforce Lorentz invariance on the theory, we 
use RI to generate the higher dimensional operators with fixed coefficients. In PQ^PT the 
fixed coefficient Lagrangian is given by 



b3^b \ + ( T»J^T U ) n 



2M 



B 



2M B 11 



-„ ib ■ s . _ 
T"— — v ■ AT,, 



M, 



v ■ A——T„ 



B 



M 



B 



(47) 



As in xPT there are additional operators with unfixed coefficients. These higher dimensional 
operators relevant for our calculation are collected in the PQ%PT Lagrangian 1 



tf Tf{A v A^'r^ k + ^(-)«7y(^(^/7j t + tf (T,T») str(A^) 



4tt/ 



+tfTf{A»A v )?T^ k + ^7yW(4)/^ + t} (T,T») tr(AM„) 

+q A Tf(v • Av A)/7% k + tl A Tf{v ■ A)f\v ■ A)/T*, k + tf (T M T") tr(« ■ Av ■ A) 



+tf (T^)tT(M+M+)+tf (T^M+T^) tr (M+)+t™ (T^)ti{M + )ti{M. 



(4? 



The coefficients of these operators tf, tf, t^ A , and tf' 1 have the same numerical values as in 
XPT. 



V. DECUPLET MASSES IN PQxPT. 

Only the masses of the octet baryons have been investigated in PQxPT ji^, HH- Here 
we calculate the masses of the decuplet baryons to NNLO in PQxPT. The chiral expansion 



1 We have omitted three PQ^PT operators of the form T M [A^, A V \T V and three of the form 
[A v , Ap\S v SfT^ since their contributions to the decuplet masses vanish to the order we are work- 
ing. Such operators identically vanish in xPT. 
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TABLE VIII: The coefficients and A^,, in PQxPT. Coefficients are listed for the decuplet 
states T, and for A^ are grouped into contributions from loop mesons with mass m^, while for 
ATj/ are grouped into contributions from pairs of quark-basis rj q mesons. 
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TABLE IX: The coefficients B^ and Bl., in PQ^PT. Coefficients are listed for the decuplet states 
T, and for 2? J are grouped into contributions from loop mesons with mass m^, while for B^, are 
grouped into contributions from pairs of quark-basis ij q mesons. 
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of the decuplet baryon masses in PQ^PT has the same form as in %PT. 

M Ti = M (//) + (/i) + Mg /2) (,x) + (ai) + • • • (49) 

The contributions at each order are similar to those in Section ITTT1 However, we must also 
include hairpin contributions from the flavor diagonal propagator, see Fig. Eland Fig. HJ 
To leading order in PQ^PT the decuplet masses are 



M^P = -7m rriT-2 o M str(m Q ), 
where the coefficients vht appear in Table HI At next-to-leading order, we have 

5H 2 



(50) 



Mr 



(3/2) 



727T/ 2 

c 2 



Y A l m l + Y A l<t>' M3 ( m ^ m 4>') 
<t> <t>4>' 



(51) 



where Tim, A, /x) is given by Eq. 



Additionally we have employed the abbreviations 



T{m^ mp, 5, /i) = [^{m^ 5, /j,), T(mp, 5, /i),f(mx, 6, fi)], 



(52) 
(53) 
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FIG. 3: In addition to the one- loop diagrams in Fig.Q] MX' 1 also receives contributions from the 
singlet (hairpins) in PQ^PT. Single and doubles line correspond to 240-baryons and 138-baryons 
respectively. The crossed dashed line denotes a hairpin propagator. The filled squares denote the 
axial coupling given in Eq. (|46jl. 

for contributions arising from the hairpin diagrams depicted in Figure H3 The sums involving 
these two functions are over pairs of flavor neutral states in the quark basis, e.g., above <p<p' 
runs over r] u r] u , rj u rj S) and rj s r] s . In this way there is no double counting. The coefficients 
and AT appear in Table IV111I and BT and BT,, appear in Table 11X1 



=■=•=■= =■ — •- 



s "X^^ ^-X- v y 'X 

/ / 
\ \ 




/X. /X. 

I I 

/ ; 

\ \ 



FIG. 4: In addition to the one-loop diagrams in Fig. H M|. also receives contributions from the 
singlet (hairpins) in PQ^PT. Single and double lines correspond to 240-baryons and 138-baryons 
respectively. The crossed dashed line denotes a hairpin insertion. Filled squares denote the axial 
coupling given in Eq. (|46|) . Filled circles denote a coupling to M.+ given in Eq. (|46|). Empty squares 
and empty circles denote insertions of fixed 1/Mb operators given in Eq. (|47j). and insertions of a 
two-baryon-two-axial coupling given in Eq. (|48[) respectively. 
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At next-to-next-to leading order, we have the contribution to the decuplet mass 



Mr 



(2) 



(Z - 1)M 



(i) 

T 



+ 4^f {l^ {m2)T + I* 2 (W)t + Str(m ^ + mrstr ( m <?) + I str ( m Q)] 2 } 
^ Cj £(m^,, /j) + ^ Cj^, C(mj, m# , y) 



+ 



27m 
(4vr/) 2 

2aAf 



(47T/) 1 



y~] C</> £(m^, /i) +y~] £(m^, , //) 



^3 ^ (tf £>J + t^J + t§ D<f) £(ni0, /i) 



(4tt/) 



+ (4^7)3 S W + *2 + *3 D w) £( m 4>i "V,/i) 



+ 



+ 



4(4tt/) 



4(4tt/) 



1— £ [(if + ^)23j + (if + tt)El + (t§ + tl A )D, 
J > 



Z(m^, M )-5m} 



H 5 



+ 



8 (4^/) 2 M B 

15 C 2 
' 16 (4tt/) 2 M s 

10?t: 2 aMstr(mQ) 
3(4vr/) 2 

3C 2 (JMStr(mQ) 

WH 2 jm 
9(4vr/) 2 

3C 2 



E4 

</> 

E B 
E A 



■7T/ N 19 4 

£(m^,,/x) + — 
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+ E^' £ ( m tf»"v>M) + Y5- A/i4 ( m *' m </'') 



Z(m*M)-^m} 



„, . 26 o 
£(m ,/i) + — m 



£(m^,m^,/i) - — M 4 (171^,171^) 



26 



<t>4>' 



^B^J (rri0, - A, m) + E ^ ("V' "V ' ~ A ' ^ 



E F * 



N 26 2 
£{1710, fx) + — 77fy 



£ {m0,m0,,ji) + —M 2 (m,p, nip) 
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2(4tt/)S 



^ ,7(m0, - A, /x) + ^ G^,, J (m^, mp , - A, /j) 



(54) 



The above expression is written quite compactly. It involves the tree-level coefficients %, 
(m 2 )T, and (mm')r, which are listed in Table HJ and the wavefunction renormalization Z, 
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TABLE X: The coefficients Cj, Cj^,, C<p, and C^r in PQ^PT. The coefficients and are 
identical for all decuplet states. The remaining coefficients are listed for the decuplet states T , 
and for Cj are grouped into contributions from loop mesons with mass m^, while for Cj^, are 
grouped into contributions from pairs of quark-basis rj q mesons. If a particular meson or pair of 
flavor neutral mesons is not listed, then the value of the coefficient is zero for all decuplet states. 
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which we find to be 
bU 2 



Z -I 



3(4vr/)2 

3C 2 
'2(4vr/)2 



£4 



C (m^, m^,/j) + ^M 2 (m^, 



5^ (m , - A, /i) + 5j 0/l 7 (m , 77V, - A - ^) 



(55) 



The new non-analytic functions arising from loop integrals in Eqs. (|54|) and ()55|) are defined 
to be 



£(m^m f ,/i) =H<f,p[C{m$,fi),£(mp,ii),£(mx,fJ>)], (56) 
£(m^,m^,//) = H (t> ^[C(m^,fi),C(m ( f >/ ,fi),C(mx,fx)}, (57) 
J (m^m^^, n) =Hw[J(m$,6,iJ$,J(mp,6,p),J(mx,S,iJL)] (58) 

and arise from the hairpin contributions shown in Figure 0] The various coefficients in the 
above sums are listed in Tables IVIIlMxTVl One may check that in the limit rrij — > m u , 
m r — > m s the xPT results are obtained from the PQxPT expressions. 

As in the xPT case, the above expressions Eqs. (joTIj) . (J53)l . and (jo^j) are quark mass 
expressions, i.e. , the meson masses in these expressions are merely replacements for the 
quark masses via Eq. (0J). To utilize these expressions for chiral extrapolations of lattice 
QCD data, one must write the quark masses in terms of the lattice meson masses. To work 
to 0(m 2 q ) means the LO contributions in Eq. (joHjl must be written in meson masses out to 
NLO. 



VI. SUMMARY 



We have calculated the masses of the decuplet baryons in the isospin limit of three-flavor 
xPT and have also derived the decuplet masses in the analogous partially quenched theory. 
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TABLE XI: The coefficients D^, D^,, D^, and in PQ%PT. The coefficients and 

are identical for all decuplet states. The remaining coefficients are listed for the decuplet states 
T, and for are grouped into contributions from loop mesons with mass m^, while for D^, are 
grouped into contributions from pairs of quark-basis rj q mesons. If a particular meson or pair of 
flavor neutral mesons is not listed, then the value of the coefficient is zero for all decuplet states. 
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TABLE XII: The coefficients and -E^/ in PQ%PT. Coefficients are listed for the decuplet states 
T, and for are grouped into contributions from loop mesons with mass m^, while for Ej^, are 
grouped into contributions from pairs of quark-basis rj q mesons. 
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We have kept all relevant terms to 0(m q ), including terms fixed by reparameterization in- 
variance that ensure the Lorentz invariance of the heavy baryon effective theory. Knowledge 
of the low-energy behavior of QCD and PQQCD is crucial to properly extrapolate lattice 
calculations from the light quark masses used on the lattice to those in nature. 

Working to 0{rri^) in the chiral expansion forces the introduction of a large number 
of low-energy constants. For any predictive power, the low-energy constants must be fit 
from experiment or lattice results. Lattice calculations will eventually allow first principles 
determination of these constants and thus predictions from QCD. In the foreseeable future, 
partially quenched simulations will enable these rigorous predictions. Our PQ%PT results 
for decuplet baryon masses are required for the proper extrapolation of PQQCD lattice data 
and hence for physical predictions from QCD. As the decuplet baryons are resonances, a 
short procedural comment is in order. For large enough pion masses, the decuplet states 
are stable to strong decays and can be studied on the lattice, cf. our expressions are all 
real valued for these pion masses. Thus for m n ~ 300 MeV, 2 where one still trusts the 
chiral expansion, the decuplet properties can be calculated on the lattice. One then uses the 
expressions we have derived to fit the low-energy constants. The expressions at next-to-next- 



2 This is a conservative estimate. Volume effects will modify the decay threshold of the decuplet. 
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TABLE XIII: The coefficients Fj and Fj^, in PQ%PT. Coefficients are listed for the decuplet states 
T, and for Fj are grouped into contributions from loop mesons with mass m^, while for Fj^ are 
grouped into contributions from pairs of quark-basis rj q mesons. 
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to-leading order are necessary for reducing uncertainty in the chiral extrapolation. Armed 
with the low-energy constants, one can then make predictions for the decuplet resonances, 
e.g. their decay widths can be found via 

3m(M T ) = -^( A2 -^) 3/2 ^ (59) 
- ^ A - ml | [ 7M m T + 3(<j m - a M )tr{m Q )\ B T „ + | G?} . (60) 

To reiterate: the low-energy constants appearing in %PT appear in PQxPT and by fitting 
them using PQQCD lattice calculations one can make QCD predictions. The decuplet 
baryon masses are no exception; our PQxPT results exhibit a smooth limit to x?T as the 
sea quark masses vary. Furthermore, lattice simulations with unphysically stable decuplet 
states can be used in conjunction with %PT and PQxPT to predict properties of the physical 
resonances. 
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TABLE XIV: The coefficients and G^, in PQxPT. Coefficients are listed for the decuplet 
states T, and for are grouped into contributions from loop mesons with mass m^, while for 
GT,, are grouped into contributions from pairs of quark-basis rj q mesons. 
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